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ABSTRACT 



Context. Rotational mixing in massive main sequence stars is predicted to monotonically increase their surface nitrogen abundance 
with time. 

Aims. We use this effect to design a method for constraining the age and the inclination angle of massive main sequence stars, given 
their observed luminosity, effective temperature, projected rotational velocity and surface nitrogen abundance. 

Methods. This method relies on stellar evolution models for different metallicities, masses and rotation rates. We use the population 
synthesis code STARMAKER to show the range of applicability of our method. 

Results. We apply this method to 79 early B-type main sequence stars near the LMC clusters NGC 2004 and Nil and the SMC 
clusters NGC 330 and NGC 346. From all stars within the sample, 17 were found to be suitable for an age analysis. For ten of them, 
which are rapidly rotating stars without a strong nitrogen enhancement, it has been previously concluded that they did not evolve as 
rotationally mixed single stars. This is confirmed by our analysis, which flags the age of these objects as highly discrepant with their 
isochrone ages. For the other seven stars, their nitrogen and isochrone ages are found to agree within error bars, what validates our 
method. Constraints on the inclination angle have been derived for the other 62 stars,with the implication that the nitrogen abundances 
of the SMC stars, for which mostly only upper limits are known, fall on average significantly below those limits. 
Conclusions. Nitrogen chronology is found to be a new useful tool for testing stellar evolution and to constrain fundamental properties 
of massive main sequence stars. A web version of this tool is provided. 
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1. Introduction 

The age determination of stars is an important objective, and sev- 
eral methods are currently being used. For star clusters, which 
presumably contain stars of similar ages, comparing the main 
sequence turn-off point in a color-magnitude or Hertzsprung- 
Russell diagram with stellar evol ution models can predict the 
age of all stars in the cluster (Hial l2002l lBuonannol ll986). For 
individual stars, ages can be derived by comparing their loca- 
tion in the HR-diagram with isochrones computed from stellar 
evolution models. In the envelope of low mass main sequence 
stars lithium is depleted, resulting in a monotonous decrease of 
the surface lithium abundance with time, which can be used to 
determine their age (Soderblom 2010; Randich 2009). 

Massive stars fuse hydrogen to helium via the CNO cycle. 
The reaction 14 N (p, y) I5 has the lowest cross section, result- 
ing in an increase of the nitrogen abundance compared to car- 
bon. The evolution of massive stars is furthermore thought to be 
affected by rotation, which may cause mixing processes includ- 
ing large scale circu lations and local hydrodynamic instabilities 
(Enda l&Sofialll978l) . The nitrogen produced in the core can be 
transported to the surface when rotational mixing is considered, 
causing a monotonic increase of the surface nitrogen abundance 
with time dHeger & Langerll2000t lMaederll2000t IPrzvbilla etaD 
l2O10h . Here, we suggest to utilize this effect to constrain the ages 
of stars. 

Based on detailed stellar evolution models, we present a 
method which reproduces the surface nitrogen abundances of 
massive main sequence stars as a function of their age, mass, 



surface rotational velocity and metallicity. Three stellar evolu- 
tion grids for the Milky Way (MW), the Large Magellanic Cloud 
(LMC) and the Small Magellanic Cloud (SMC) are used, pro- 
viding stellar models of different metallicities, masses and ro- 
tational velocities and include the effects of rotational mixing 
(iBrott et al.ll20T Tab . The mixing processes in these stellar evo- 
lution models were calibrated to reproduce the massive main 
sequence stars observed within the V LT-FLAMES Survey of 
Massive Stars dEvans et al.ll20Q5L 12006b . 

The abundance analysis of massive stars is often focused on 
stars with small v sin i va lues, as this allows t o obtain abundances 
with higher accuracies (Przvbilla et al.l 120101) . However, those 
objects are not well suited for nitrogen chronology, as they are 
either slow rotators, or have a low inclination angle. In the lat- 
ter case, the age can not be well constrained, since a large range 
of true rotational velocities is possible. Many stars analyzed in 
the VLT-FLAMES Survey, however, are particularly suited as 
targets for nitrogen chronology, since this survey did not con- 
centrate on the apparent slow rotators. 

On the other hand, results from the VLT-FLAMES Survey 
identified two groups of early B-type stars which could not be 
reconciled by the evolutionary models of single stars with rota- 
tional mixing dHunter et al.ll2008at IBrott et al.ll201 lbh . The first 
group consists of slowly rotating stars which are nitrogen rich, 
to which the method of nitrogen chronology is not applicable. 
The second group contains relatively fast rotating stars with only 
a modest surface nitrogen enhancement. While suc h stars are 
in prin ciple predicted by models of rotating stars (Brot t et al.l 
1201 lal) . they were not expected in be found in the VLT-FLAMES 
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Survey due to the imposed V-band magnitude cut-off, since un- 
evolved main sequence stars are visually faint, and the visually 
brighter more evolv ed rapid rotators a r e expected to be stro ngly 
nitrogen enriched (Brot t et al.1l201 lbl) . iBrott et all (1201 lbb sug- 
gested that either these stars are the product of close binary evo- 
lution, or that rotational mixing is much less efficient than pre- 
dicted by the stellar evolution models. 

The nitrogen chronology is able to shed new light on 
this issue. In cont rast to the population synthesis pursued by 
dBrott et al.ll201 lbl) . it does not use statistical arguments but is 
applied to each star individually. Its results can therefore be eas- 
ily compared with classical age analyses, e.g., using isochrone 
fitting in the HR diagram. In this paper, we apply the nitrogen 
chronology to all LMC and SMC stars from the VLT-FLAMES 
Survey which have a projected rotational velocity larger than 
lOOkm/s. 

We organize our paper as follows. In Sect. |2]our method of 
nitrogen chronology is presented. Sect. |3]describes our analysis 
of the VLT-FLAMES sample. A discussion and evaluation of our 
new method is given in Sect.|4]and|5] Sect.|6]will round off with 
our conclusions. 

2. Method 

The enrichment of nitrogen at the surface of a massive star dur- 
ing its evolution on the main sequence depends on its mass and 
rotation rate. In Fig. [TJthe surface nitrogen abundance is plotted 
against the fractional main sequence lifetime t = f/rn for dif- 
ferent initial masses and surface rotational velocities. Here, th 
is the hydrogen burning lifetime of a stellar evolution sequence. 
A parameterization of th as a function of mass and surface ro- 
tational velocity, for all three metallicities considered here, is 
given in Appendix IA.2I For our method, we define "zero-age" 
as the time when the mass fraction of hydrogen has decreased 
by 2% in the center due to hydrogen core burning. This ensures 
that the stellar models are still unevolved but already thermally 
relaxed. The surface nitrogen abundance N is given by 

Ar = log(|) + 12, (1) 

with the surface number density of nitrogen and hydrogen atoms 
N and H, respectively. 

Figure [TJ shows that the nitrogen abundance rises smoothly 
over time starting with the surface nitrogen abundance at the 
ZAMS depending on the metallicity. It reaches its maximum 
value at the end of the main sequence evolution. The right panel 
of Fig. [TJ depicts the mass dependence of the nitrogen enrich- 
ment at the surface. For a given initial surface rotational ve- 
locity, the models show a faster enrichment at the surface for 
higher initial mass. Furthermore, the faster the initial surface ro- 
tational velocity for a given mass the greater is the effect of rota- 
tional mixing which leads to faster enrichment and a higher sat- 
uration value dHeeer & Laneerl 120001 iMevnet & Maederll2000t 
Maed er & Mevnetll2001l)~ 

To constrain the age of an observed massive main sequence 
star using its surface nitrogen abundance, we derived an equation 
which fits the surface nitrogen abundance for a given mass, sur- 
face rotational velocity and metallicity to the models of rotation- 
ally mixed massive stars. In order to reduce interpolation errors 
in constructing this equation, stellar evolution grids are required 
which have a dense spacing in the grid parameters initial mass M 
and initial equatorial rotational velocity Vzams- Presently, stellar 
evolution grids ava ilable for this purpose are those presented in 
iBrottet al.ld2011ah . 



The surface nitrogen abundance monotonically increases 
over the fraction of the main sequence lifetime (cf. Fig. [TJ. This 
can be reproduced well with an ansatz as 

N - a + c ■ y d (b ■ t) . (2) 

The parameters a, b, c, and d are functions of the initial surface 
rotational velocity Vzams an d the initial mass M of a star. The in- 
complete gamma function jd (b ■ t) depends on the product b ■ t, 
with the fraction of the main sequence lifetime r and the param- 
eters b and d. For calculating the surface nitrogen abundance ac- 
cording to Eq. © the incomplete gamma function is normalized 
by the complete gamma function T as 

yd{b-T) = -±— f v' ; l e \U. (3) 

r(d) Jo 

Considering the properties of the incomplete gamma func- 
tion it is possible to make some statements about the parameters 
a, b, c, and d. The values of yd range from to 1, so c is re- 
sponsible for the height of the curve. The parameter a describes 
the initial surface nitrogen abundance, which is independent of 
mass and surface rotational velocity. The parameter b describes 
the speed of the nitrogen enrichment i.e. the steepness of the 
curves in Fig.Q] 

The parameters b and c depend on initial mass and rotational 
velocity. By fitting Eq. (f2| to the data of the stellar evolution 
models, both parameters can be determined for different values 
of M and vzams- Considering the parameter d, the best fits of 
Eq. (f2]i to the data of the three stellar evolution grids are obtained 
by a linear dependence of d on the surface rotational velocity 
vzams (see Appendix IA.U . 

Nitrogen profiles calculated by using Eq. (f2]i can be useful 
for binary population studies, where a quick knowledge of the 
surface nitrogen abundance is required. 

2.1. Accuracy 



Table 1. Parameters of the example sequence discussed in the 
text and displayed in Fig. [2] 



parameter 


Model 2 


Model 1 


M[M Q ] 


30 


35 


Vzams [km/s] 


169 


168 


Z 


Zlmc 


Zlmc 


a 


6.9 


6.9 


b 


12.8 


14.6 


c 


0.6 


0.6 


d 


6.0 


6.0 



Two stellar models are chosen to indicate the quality of the 
presented method to calculate the surface nitrogen abundance. 
All parameters of Eq. (f2) are calculated from Eq. (IA.1I - IA.4> 
given in Appendix lA.lK TablefTt. Using these values, the surface 
nitrogen abundance of the two sequences as a function of r is 
calculated and shown as blue solid lines in Fig. [2] It is compared 
to the nitrogen abundance evoluti on which results f rom the de- 
tailed stellar evolution models of Brott et al ] (l2011al) . 

All calculated surface nitrogen abundance functions fit the 
detailed models very well for low values of r, but show differ- 
ent levels of deviations towards the end of the main sequence 
evolution. In Fig.|2l a typical fit is shown on the left panel. At all 
times, the fit result deviates by less than 0. 1 dex from the detailed 
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Fig. 1. For different models of the stellar evolution grid of lBrott et al.l(l20 1 lab using LMC metallicity, the surface nitrogen abundance 
is plotted against the fractional main sequence lifetime. The left panel shows the influence of the initial rotational velocity vzams 
for a given mass of M — 25 M Q . The right panel depicts the evolution of the surface nitrogen abundance for several masses at 
vzams ~ 280 km/s. 
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Fig. 2. Surface nitrogen abundance as functi on of the fractiona l main sequence lifetime for two specified stellar models (see legend) 
from the detailed stellar evolution grids of lBrott et al.l (120 1 lab (red lines), compared to the results of our fit (blue solid lines). The 
blue dashed lines represent error ranges of +0.1 and +0.2 dex for the nitrogen abundance. While the left panel shows a typical fit, 
the right panel depicts our worst case. 



model prediction. The right panel of Fig.|2]shows our worst case. 
Even in this case, the fitted surface nitrogen abundance deviates 
by more than about 0. 1 dex from the value of the detailed stellar 
evolution model for a fraction of the time. 



Numerical noise in the detailed stellar evolution models does 
in rare cases lead to a non-monotonic behaviour of the nitrogen 
evolution as function of the initial rotation rate or as function of 
mass (cf., left panel of Fig. 1). These features are smoothed out 
by our analytic fitting method, which guarantees monotonic be- 
haviour in the considered parameter space (Sect. 2.2). Therefore, 
a deviation of our nitrogen prediction from that of the detailed 
stellar evolution models does not necessarily flag a problem 
of our method. In summary, we consider the error in the pre- 
dicted nitrogen abundance introduced by our fitting method to 
be smaller than 0. 1 dex. In the present paper, the surface nitro- 
gen abunda nces of the analyze d stars have an observational error 
of ±0.2 dex dHunter et alJl2007h . 



2.2. Validity range 

As the lowest mass of the underlying stellar evolution models 
used in our fitting method is 5 M Q , this defines the lower mass 
limit for its applicability. Towards very high masses, effects of 
stellar wind mass loss produce features which are not captured 
by our method above an upper limiting mass, which depends on 
metallicity (see Tabled. 

The stellar evolution models of Br ott et al. for very 

rapidly rotating stars may undergo so called quasi-chemically 
homogeneous evolution. The behaviour of those models is also 
not captured by our fitting method. This problem is avoided 
when only initial rotational velocities below 350 km/s are con- 
sidered, which thus defines the upper velocity limit for our 
method. This also ensures that 2D effects of very rapid rotation 
on the stellar atmosphere do not play a role here. 

As outlined in Sect. IA.21 in applying our method we assume 
that the stellar surface rotation velocity remains roughly constant 
during the main sequence evolution. This is the predic tion of the 
underlying stellar evolution models (iBrott et al.ll201 lab for stars 
with initial masses below 30 M Q . For more massive stars, this 
feature still holds as long as the stellar effective temperature is 
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Table 2. Range of initial stellar parameters for which our method 
is validated, for the three considered metallicities. In addition, 
for stars more massive than 30 M , the effective temperature 
should be higher than 25 000 K. 



Galaxy 


M 


fZAMS 




M Q 


km/s 


MW 


5-35 


0-350 


LMC 


5-45 


0-350 


SMC 


5-50 


0-350 



larger than 25 000 K. For cooler temperatures, mass loss leads to 
a noticeable spin down of the stars. Our method is therefore valid 
in the ranges given in Table [2] with a restriction that for stars 
more massive than 30 M , their effective temperature should be 
higher than 25 000 K. 

3. Application 

3.1. Stellar sample 

The obtained equations to calculate the surface nitrogen abun- 
dance for given initial mass and surface rotational velocity are 
applied here to constrain the properties of observed stars under 
the hypothesis that the underlying stellar models describe their 
evolution correctly. The V LT-FLAMES Survey of Massive Stars 
dEvan s et al. 2005, 2006) has been chosen to provide the obser- 
vational data. 

The VLT-FLAMES Survey is a spectroscopic survey of hot 
massive stars performed at the Very Large Telescope of the 
European Southern Observatory. Over 700 O and early B type 
st ars in the Gal a xy and the Magellanic Clouds w ere observed. 
In IHunter et alj d2008bl) (Magellanic Clouds) and iDufton et al] 
(I2006I) (Galaxy) the spectral data has been eval uated using the 
non-L TE TLUSTY model atmosphere code dHubenv & Land 
1 1995b to gain atmospheric parameters and projected rotational 
velocities of the observed stars. The surface nitroge n abundances 
were derived using the nitrogen absorption lines (iHunter et alJ 
120091). The evo l utionary masses of the stars were determined 



Hunte r et al. (2008b) using the stella r evolu tion models of 



iMevnet et all Id 1994 an d Schaer er et all dl993l) for LMC ob- 
jects and in iMevnet et alj d!994l) and ICharbonnel et al.l d 1993b 
for SMC objects. The derived evolutionary masses are in good 
agreement with evolutionary masses determined from the stellar 
evolution models used here. 

For our analysis, we selected stars from the FLAMES LMC 
and SMC samples by several conditions. It is necessary to 
know the surface nitrogen abundance, the evolutionary mass and 
the (projected) surface rotational velocity of a star to use our 
method. We disregarded all sample stars for which this is not the 
case, or for which one of the parameters falls outside the validity 
range of our method (Sect. 12.2) . We further discarded all stars 
with a projected surface rotational velocity below 100 km/s, as 
most of them are expected to be intrinsically slow rotators for 
which no significant age constraints from our method could be 
expected. In this way, 44 LMC and 35 SMC stars remain on the 
list, which we call "our" sample in the following (Tables lBTI and 
lB~2l . 

While we need to exclude stars with rotational velocities 
above 350 km/s, it is possible that stars which rotate faster than 
this are present in our sample, since we only know the projected 
rotational velocities of our sample stars. The largest v sin(;) of 
all stars in our sample is 323 km/s. A p riori, given the dis tribu- 
tion of rotational velocities presented in dBrott et al .11201 lbh (see 



Table© of stars in the FLAMES LMC and SMC samples as well 
as a uniform distribution of orientations, 4.8% of all stars with 
vsin(/) < 350 km/s have surface rotational velocities v rot higher 
than 350 km/s. Such fast rotating stars rapidly enrich nitrogen at 
the surface. Within the first 10% of their main sequence lifetime, 
their surface nitrogen abundance increases above N — 8 dex, 
in the LMC. For a rough estimate, assume that approximately 
4.8% of all stars with vsin(/) < 350 km/s could have a surface 
rotational velocity above 350 km/s with 10% of their main se- 
quence lifetime corresponding to a surface nitrogen abundances 
N < 8 dex. As we have no stars with N > 8 dex in our sample, 
0.48% of all stars within our sample could have a surface rota- 
tional velocity above 350 km/s, which we consider negligible. 

The largest v sin(/) of our stars corresponds to about 50% of 
break-up velocity. Consequently, the largest centrifugal acceler- 
ation in the line of sight amounts to one quarter of the gravita- 
tional acceleration. In other words, the measured log g will be 
decreased by at most about 0.12 dex, and on average by less 
than half that value. This is smaller than the intrinsic measure- 
ment error of logg of about +0.15 dex. It is also small com- 
pared to the range of gravities of the sample stars, which is 
4.4 dex < logg < 3.2 dex. For this reason, we can neglect the 
correction of the gravity due to rotation. 

3.2. Defining Class 1 stars 

To constrain the age of an observed star through our method, its 
evolutionary mass, projected surface rotational velocity, metal- 
licity and the surface nitrogen abundance need to be known. 
We want to mention, that in the case of diffe rences betwee n 
evolutionary and spectroscopic masses (Moki em et al.l [2006). 
this method will lead to different results when using the spec- 
troscopic mass for deriving the nitrogen profile. The effect of 
changes in mass during the main sequence evolution, for stars 
within the validity range given in Table|2] on the calculated sur- 
face nitrogen abundance profile is smaller than the typical obser- 
vational error of ±0.2 dex. The mass, derived from the observa- 
tional data is therefore used as the input value for M in Eq. (0. 

If the true rotational velocity of the star is known (which is 
rarely the case, and not so for any of our sample stars), it can 
be used as input into Eq. (0, assuming that it still represents 
its initial rotational velocity. For a given metallicity the surface 
nitrogen abundance as a function of time is then calculated us- 
ing the equations in Appendix IA.1I Next, the fractional main 
sequence lifetime is determined by the condition that the calcu- 
lated surface nitrogen abundance matches the observed one. By 
considering the main sequence lifetime of stars as function of 
mass and surface rotational velocity (cf. Appendix lA.2K the age 
of the star is obtained. 

If the true rotational velocity of a star is not known, but only 
its projected value v sin(/), we can compute the nitrogen surface 
abundance as function of t for a model with the appropriate mass 
and metallicity, and with a rotational velocity equal to the pro- 
jected rotational velocity of our target star. Since the projected 
rotational velocity is smaller than the true velocity, the nitrogen 
enrichment will be slower than in a model rotating with the true 
rotational velocity. We therefore obtain lower limits to the ex- 
pected nitrogen abundance, for all times. The result of the anal- 
ysis depends now on the measured nitrogen abundance of our 
star. If it is lower than the largest computed lower limit (which 
is obtained for r = 1), we can constrain the upper limit to the 
current age of the star. 

We define stars of Class 1 by the condition that their ob- 
served surface nitrogen abundance, or its observationally deter- 
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mined upper limit, is lower than the maximum (final) calculated 
surface nitrogen abundance according to Eq. (0, when adopting 
the observed v sin(/) as rotational velocity in Eq. (O. 



Table 3. Parameters used in our population synthesis calculation 
(cf.Fig.g). 



X 



7.8 



7.6 



7.4 



7.2 



1 1 1 1 

Class 1 v zams= 267 km/s - 


— -"""^ i i i 


N11-95 

vsin(i)=267 km/s 
log(N/H)+1 2=7.46 



0.2 0.4 0.6 0.8 1 

fraction of the MS lifetime x 

Fig. 3. The analysis of the Class 1 stars N 11-95 (LMC), which 
has an observed vsin(/) of 267 km/s. The blue solid line shows 
the calculated surface nitrogen abundance as function of the frac- 
tional main sequence life time t of an LMC star with the mass of 
N 1 1-95 and assuming a rotational velocity of 267 km/s. The red 
dotted line indicates the measured surface nitrogen abundance of 
N 1 1-95. The black vertical arrow gives the derived upper limit 
on t, i.e. r < 0.43. 

Figure [3] illustrates our method at the example of the star 
N 11-95 from our LMC sample, with M = 15 M e , vsin(/) = 
267 km/s and N = 7.46 dex. The observed abundance of ni- 
trogen at the surface equals the calculated nitrogen abundance 
(which is a function of t) at a value of r = 0.43. Considering a 
main sequence lifetime of 11.4 Myr, this leads to an upper age 
limit of 4.9 Myr for this star. 

For some of the sample stars we only know upper limits to 
the surface nitrogen abundance. Since the nitrogen upper limit 
is larger than the true nitrogen abundance of the sample star, 
using the upper limit to the nitrogen abundance leads to an upper 
age limit of the observed star in the same way as for stars with 
determined nitrogen abundances. 

3.3. Population synthesis of Class 1 stars 

It is interesting to consider the location of Class 1 stars in the 
so called Hunter-diagram, where the surface nitrogen abundance 
of stars is plotted versus their projected rotational velocity. Since 
according to the stellar evolution models fast rotators are achiev- 
ing higher surface nitrogen abundances, stars are likely to fall 
into Class 1 when they have a large v sin(7) (which makes their 
maximum nitrogen abundance large) as long as their measured 
nitrogen abundance is not too high. Conversely, nitrogen rich ap- 
parent slow rotators will not fall into Class 1 . 

To study this, we perform a population synthesis simulation 
for rotating single massive main seque nce stars of LMC metal- 
licity using the code STARMAKER ( Brott et alJfeOllbh . with 
parameters as summarized in Table [3] 

To verify the Class 1 character of a simulated star with a 
given age, its mass and v sin(/) value have been used to compute 
its surface nitrogen abundance at the end of its main sequence 
evolution according to Eq. (0. If this value is higher than its 



parameter 


star formation rate 


constant 


considered time interval 


35 Myr 


velocity distribution 


Gaussian distribution (Brott et al. 201 lb) 




(<x= 141 km/s, yu=100 km/s) 


velocity range 


0-570 km/s 


mass distribution 


Safpeter ("1955) IMF 


mass range 


5 - 30 M Q 


stellar model grid 


LMC (Brott etal. 2011a) 



actual surface nitrogen abundance, it is taken to be of Class 1 . As 
all stars in the LMC and SMC sample have masses below 30 M G 
the maximum mass considered in this simulation was chosen as 
30 Mq. Only stars with vsin(z) > 100 km/s are evaluated. 



+ 




b 10000 



1000 



100 



10 



150 200 250 300 350 400 
vsin(i) [km/s] 

Fig. 4. Hunter-diagram of Class 1 stars of a synthetic population 
of main sequence s tars obtained using the STARMAKER code 
(Brot t et aT[|201 lbh . The surface nitrogen abundance is plotted 
versus the projected surface rotational velocity. The colors indi- 
cate the number of Class 1 stars in a given pixel. The absolute 
scali ng is arbitrary. The Hunter-diagram is divided into 5 areas 
as in (Brott et al. 20 1 1 hi) to simplify the later analysis. 

Figure |4] shows the number density of simulated Class 1 
stars in the Hunter-diagram. Stars are counted within intervals 
of Av sin(/) = 10 km/s and AN = 0.05 dex. The number of mod- 
els in the intervals are indicated by the color coding. 

Figure [4] illustrates for which stars the nitrogen chronology 
method can be used to derive an upper age limit. As expected, 
preferentially stars in the lower right corner of the Hunter- 
diagram can be evaluated as Class 1. For stars in the upper 
left corner of the Hunter-diagram we will be able to constrain 
their inclination angles (see below). Finally, we note that if the 
inclination angle of a star is observationally determined (e.g. 
through a pulsation analysis, or a binary membership), its ni- 
trogen age can be determined regardless of its location in the 
Hunter-diagram. 

The majority of Class 1 stars in our sample are located 
in the LMC, only one star belongs to the SMC sample (see 
Appendix [B]). The reason is that nitrogen abundances of SMC 
stars are lower than those of comparable LMC (or Galactic) stars 
dHunter et al.l l2009). such that their abundance can not be deter- 
mined if it is below a certain threshold which increases with the 
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projected rotation velocity. That means, Box 1 in the Hunter- 
diagram of t he SMC early B sta r sample is essentially empty 
(See Fig. 7 in lHunter et ail d2009l) ). 

Furthermore, a comparison of Fig. |4] with Fig. 9 of 
iBrott et ail (1201 lbl) shows that only very few stars from the con- 
sidered sample are expected in Box 1, due to the visual magni- 
tude cut-of f impos ed in the sample selection. For this reason, 
Brott et al.l (1201 lbl) concluded that those stars of this sample 
which do appear in Box 1 are not evolving as rotationally mixed 
single stars. 

3.4. Defining Class 2 stars 

We define Class 2 stars by the condition that the observed ni- 
trogen abundance, or its observationally derived upper limit, 
is higher than the maximum (final) calculated surface nitrogen 
abundance according to Eq. (f2]) when adopting the observed 
v sin(z') as rotational velocity in Eq. ([2]). Stars are thus either of 
Class 1 or of Class 2. For Class 2 stars, the derived upper limit 
to the current age is the main sequence lifetime, which yields no 
age constraint. 

However, for a Class 2 stars with a determined nitrogen 
abundance, we can find the minimum required rotational veloc- 
ity Vmin to achieve this nitrogen value for the given mass of the 
star. Since the surface nitrogen abundance is monotonically in- 
creasing with time, the observed nitrogen abundance will then 
be produced by a model star at the end of the main sequence 
evolution (t = 1), rotating with this minimum velocity. By using 
the minimum rotational velocity, the maximum inclination an- 
gle which is compatible with the observed nitrogen abundance, 
which we designate as nitrogen inclination z'n, can be derived for 
the star. The sine of z'n can be calculated as 



sin (z' N ) = 



v sin (z) 



(4) 
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Fig. 5. The determination of the nitrogen inclination angle of a 
Class 2 star is shown for the example NGC 2004-30. The red 
dotted line indicates the observed surface nitrogen abundance 
and the blue solid line represents the calculated surface nitrogen 
abundance as a function of the fraction of the main sequence 
lifetime using vzams = 123 km/s in Eq. (0. The blue dashed 
line is calculated with a value of the surface rotation velocity of 
vzams = 173 km/s. This corresponds to an inclination angle of 
sin(z'N) = 0.71. 



Figure [5] illustrates the analysis of Class 2 stars at the exam- 
ple of the LMC star NGC 2004-30. To calculate the evolution 
of the surface nitrogen abundance, M = 19M and vzams = 
123 km/s have been used. The observed nitrogen value of N — 
7.43 dex is higher then the predicted one during the whole main 
sequence evolution. The observed nitrogen value can be reached 
with a minimum rotational velocity of v m ; n = 172 km/s, which 
leads to sin(z'N) = 0.71 (blue dotted line), i.e. z'n = 45.6°. 

For Class 2 stars where only an upper limit to the surface 
nitrogen abundance has been determined, a nitrogen inclination 
can be determined in the same way, by adopting the upper limit 
as nitrogen abundance. However, in this case, the nitrogen incli- 
nation does not provide an upper limit to the true stellar incli- 
nation, since lower nitrogen surface abundances would lead to 
larger nitrogen inclinations. Their values are nevertheless inter- 
esting, as shown in Sect. 5 for the Class 2 stars in the SMC. 

4. Results: Class 1 stars 

4. 1 . Isochrone ages 

To evaluate the age constraints obtained for Class 1 stars via the 
nitrogen abundances they are compared to the results from the 
well proven method of estimating the age of a star through its 
position in the Hertzsprung-Russel (HR) diagram. The idea is 
to compare effective temperature and luminosity of the star with 
isochrones plotted in the HR diagram. 

We computed isochron es using non-rot ating models of the 
stellar evolution grids of (IBrott et al.l 12011 at up to an age of 
50 Myr with a step size of 0.2 Myr. Comparing the isochrones 
with the position of a star the best fitting isochrone defines its 
isochrone age fi SO chrone as shown exemplary for the star N 1 1-95 
in Fig. [6] The observational error bars imply an uncertainty in 
the isochrone ages of about + 1 Myr. 
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Fig. 6. Isochrone fitting in the HR diagram for the LMC star N 
11-95, which is indicated by the red point with error bars. The 
red, green and blue lines are isochrones from non rotating stel- 
lar models for different ages. The green line defines the most 
likely age of the star while red and blue represent maximum and 
minimum possible current age when considering uncertainties in 
measured effective temperature and luminosity of the star. 

All stars of Class 1 have been analyzed this way to obtain 
their isochrone age. Since non-rotating models were used, rota- 
tional mixing is not considered in estimating the isochrone ages 
of the stars. However it leads only to negligible errors in the 
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mass and rotational velocity range considered here (see Fig. 
in Appendix IA. lb . All isochrone ages are given in Tables IBTTI 
andlfel 

Assuming that the isochrone age reproduces the current age 
of a star well, our method allows to derive its inclination an- 
gle. Due to uncertainties in the isochrone age fi SO chrone, and when 
using Eq. ( IA.5b to derive the main sequence lifetime th, it is 
possible that r H < fisochrone. In this case we limit fi SO chrone = t h 
(occurred in 3 out of 16 cases). To derive the inclination angle, 
the isochrone age needs to be smaller than the upper limit to the 
current age using the nitrogen abundance. Amongst the stars of 
Class 1 (the method is similar to deriving the inclination angle 
for Class 2 stars, see Sect. 13.4b . a few stars fulfill this condition. 
The derived inclination angles are given in Tables I5TI and IB. 21 
In particular, we find that the most likely inclination angles of 
stars L15 (N 1 1-120) and L42 (NGC 2004-107) are 66° and 62°, 
respectively, while L6 (N 11-89) and L8 (N 11-102) are seen 
nearly equator-on. 

4.2. Nitrogen ages versus isochrone ages 

The results of the age constraints through nitrogen are sum- 
marized in Tables |B~T1 and lB~2l for all considered stars. Fig. [7j 
shows the upper age limits of our sample stars obtained via ni- 
trogen chronology ?n, plotted against the ages obtained through 
isochrone fitting, fi SO chrone- Only stars of Class 1 are depicted. 
Due to the distinction between Class 1 and Class 2, the stars 
in Fig. [7] are those which have preferentially low surface nitro- 
gen abundances and high projected surface rotational velocities. 
Figure [9] shows the distribution of Class 1 stars in the HR dia- 
gram. 

Since the projected rotational velocity provides a lower limit 
to the true rotational velocity of a star, the surface nitrogen abun- 
dance results in an upper age limit, i.e., the star is younger or as 
old as the obtained nitrogen age. Therefore, the upper age limit 
derived from the surface nitrogen abundance should be larger 
than the age determined using isochrones, under the assumption 
that the age of a star is approximated well by the second method. 
This means that stars should be placed below the green dotted 
line in Fig. [7] which corresponds to fi soc hrone = t$. 

Considering that the error in mass is of the order of 10-25% 
(Mok iem et al.ll2006f) — which is insignificant given the weak 
mass dependence of the nitrogen enrichment in the considered 
mass range (F ig.[Tji — and the error in v sin(/) is only about 10% 
dHunter et al.l l2008b). th e error of ± 0.2 dex for the surface nitro- 
gen abundance ( Hunter et al.ll2009h is indeed the dominant one. 
To estimate its influence, we repeat the analysis for each star, 
comparing the surface nitrogen abundance according to Eq. d2J 
with an abundance which is 0.2 dex higher or lower than the 
measured surface nitrogen abundance. The lower nitrogen abun- 
dance leads to a decrease while the upper value results in an 
increase of the upper age limit for the star. The uncertainties of 
the upper age limit due to the uncertainties in the surface nitro- 
gen abundance are given in Tables IBTTI and IB~2l and are plotted 
as horizontal error bars in Fig. [TJ 

The majority of our sample stars are placed above the green 
dotted line in Fig. |7l even when considering the error in their 
measured surface nitrogen abundance. This contradicts the ex- 
pected behaviour. For example, the star N 1 1-95 (star L7 in our 
sample, see Table IB. 11 analyzed in Fig.[3]has an upper age limit 
from nitrogen chronology of 5.8 Myr. However, its isochrone 
age is 10.4 Myr. This means that, according to its position in the 
HR diagram, this rapidly rotating star (vsin(/) = 267km/s) is 



seen near its TAMS position in the HR diagram, but it has only 
a modest surface nitrogen enhancement. 

Only two out of 17 stars are found in the predicted area (N 
11-120=L15, and NGC 2004-107=L42). Considering the uncer- 
tainty of +0.2 dex in the measured surface nitrogen abundance, 
we see in Fig.|7jthree more stars (N 1 1-102=L8, N 1 1-1 16=L12, 
and NGC 330-5 1=S25) with an error bar reaching into the al- 
lowed area, and two further stars (N 1 1-89=L6 and NGC 2004- 
88=L35) with an error bar getting very close to the allowed area. 
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Fig. 7. The upper age limits obtained from nitrogen chronol- 
ogy for Class 1 stars are plotted against the ages derived using 
isochrone fitting in the HR diagram, isochrone- Stars located in 
the LMC are shown with red symbols while the SMC star is de- 
picted in blue. Stars for which only upper limits to the surface 
nitrogen abundances could be determined are marked by filled 
symbols with arrows indicating the direction of change in the 
results when the surface nitrogen abundance is smaller then the 
value of the upper limit. Circles mark single stars and triangles 
indicate probable binaries. The green dotted line is defined by 

— ^isochrone- 



4.3. Discussion 

As shown above, the nitrogen chronology clearly fails for 10 out 
of the 17 Class 1 stars in our sample. It is interesting to consider 
the location of these stars in the Hunter-diagram (Fig. [8). It turns 
out that all 10 stars are located in Box 1. In contrast, amongst 
the 7 stars for which the nitrogen chronology might work, the 
six LMC stars are found in Box 3 (L8, L15), in Box 5 (L6, L42), 
or in Box 1 but very near Box 5 (L12, L35). The SMC star S25 
can not be compared in the same Hunter-diagram. For the LMC 
Class 1 stars, also the reverse statement holds: Nitrogen chronol- 
ogy might work for all stars which are not in Box 1, or which 
are in Box 1 very near Box 5. In other words, a slight shift of the 
dividing line between Box 5 and Box 1 would result in separat- 
ing the Class 1 stars for which nitrogen chronology clearly fails 
(Box 1) from those for which it might work (outside of Box 1). 

The isochrone age of the 10 LMC stars in Box 1 for which 
nitrogen chronology fails is larger than the upper age limit de- 
rived from their surface nitrogen abundance, given their rather 
high projected rotational velocities. Due to the simplicity of our 
method the revealed conflict can only be due to the stars in Box 1 
behaving differently than predicted by the evolutiona ry models 
for rotationally mixed single stars of lBrott et al.l(l201 lab . 
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Fig. 8. Hunter-diagram of all stars shown in Fig. |7] The surface 
nitrogen abundance is plotted as a function of the projected sur- 
face rotational velocity for stars located in the LMC. Upper lim- 
its to the surface nitrogen abundance are marked by filled sym- 
bols. In addition there is a differentiation between single stars 
(circles) and probable b inaries (triangles). A subdivision into 
five boxes as defined in iBrott et a"D d201 lbl) has been done to 
simplify the analysis. The green dotted line represents the initial 
surface nitrogen abundance for LMC B stars. 
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Fig. 9. The Hertzsprung-Russel diagram is shown for several 
evolutionary models (LMC) with initial masses from 9 to 30 M 
and surface rotational velocities of approximately 200 km/s. 
Additionally, our sample with LMC-metallicity is depicted. 
Stars of Class 1 are indicated by circles and stars of Class 2 are 
represented by triangles. With respect to Fig. [8] stars in Box 1 
are shown in red, star in Box 3 in blue and stars in Box 5 in 
black. Grey dots depict the remaining not analyzed stars of the 
V LT-FLAMES Surve y (LMC) after the selection criteria applied 
in lBrottet all (1201 lbl) . 



The detailed and tailored population synthesis of the LMC 
early B type main sequence stars from the FLAMES Survey by 
Brott et all (1201 lbl) already pointed out that the stars in Box 1 
could, in a statistical sense, not be reproduced through the stel- 
lar models of single stars including rotational mixing. Our ni- 
trogen chronology analysis is not a statistical method, but rather 
considers each star individually. While it uses the sa me under- 
lying stellar evolution models as IBrott et al.l (1201 lbl) . the latter 
did not consider age constraints. In this respect, our chronology 
method provides a new test of the evolutionary models provided 



by Brott et al. (1201 lal) which is independent of the results of 
Brot t et al.1 (1201 lbl). Nevertheless, our result is very much in line 
with that of lBrott et al.l d201 lbl) . i.e., that Box 1 stars in our sam- 
ple are generally not consistent with the models of Brott et al. 
(1201 Id) . 

While the errors due to inaccuracies of our fitting method are 
mostly insignificant compared to the observational errors in ni- 
trogen, the main theoretical error source could be the systematic 
error introduced by the physics included in the underlying stel- 
lar evolution models. However, the mixing physics (convective 
overshooting and rotational mixing) in these models was specif- 
ically calibrated using the LMC sample of the FLAMES Survey. 

We are left with the conclusion that either rotational mix- 
ing does not operate in massive stars, or that the stars in Box 1 
of the Hunter-diagram (Fig. [8} are affected by physics which is 
not included in the stellar evolution models (e.g., binary effects). 
In this context, we point out that while rotating stars may well 
evolve into Box 1 in general (cf. Fig.|4]i, most of them would not 
be detected in the FLAMES survey, due to the employed mag- 
nitude cut-off (cf. Figs. 6 and 10 of lBrottet ail (1201 lbl) ). The 
origin of these stars is therefore in question. 

For example, it has been shown that binary evolution can 
place stars into Box 1 in the Hunter-diagram (Langer et al. 
2008; Ide Mink etal] 1201 lb for which the age constraints de- 
rived from nitrogen chronology do not apply. However, while 
four Class 1 stars in our sample are likely binaries (Fig. [7}, 
the remaining twelve show no indication of binarity, despite 
the fact that multi epoch spectra are available for them. The 
FLAMES Survey observing campaign was not tailo red for bi- 
nary d etection. Applying the method developed by San a et al] 
(120091) the limited number of epochs should roughly result in a 
detection probability of at least 80% for systems with periods 
up to 100 days, dropping rapidly to a detection rate of less than 
20% for periods longer than a year. Furthermore, while some 
binaries amongst the Box 1 stars may remain undetected, it is 
the pre-interaction b inaries which are more likely to be detected 
(Ide Mink et aljl20~Pll) . With the data at hand, it therefore seems 
unlikely that we can evaluate the possibility that binarity by itself 
can completely explain the situation. 

When we apply nitrogen chronology to LMC stars outside of 
Box 1 in Fig. [8] — i.e., to the only stars whic h can be expected to 
be compatible with the single star models o f IBrott etal.1 (120 Hal) 
— the result is positive. While there are only four stars (L6, L8, 
L15, L42), their nitrogen age is consistent with their isochrone 
age: They are either in the allowed part of Fig. |7](L42, L15), or 
nearest to its borderline (L6, L8). As explained above, a slight 
shift in the borderline between Boxes 1 and 5 in Fig.|8]includes 
two more stars (L12, L35) into this consideration. Together with 
the SMC star S25, for which can not be considered in Fig. [8] 
there are seven stars for which the nitrogen chronology yields 
promising and consistent results. 

In summary, we conclude that our new method give results 
for Class 1 stars which are consistent with those from classical 
isochrones for all stars in the subgroup where this could be ex- 
pected. Furthermore, it flags contradic tions for al l stars i n the 
second subgroup which, according to IBrott et aTl (1201 lbl) . are 
not thought to have evolved as rotationally mixed single stars. 
For stars of both subgroups, the application of nitrogen chronol- 
ogy is fruitful and provides insights which are not obtainable 
otherwise. 

Obviously, it would be desirable to apply nitrogen chronol- 
ogy to a stellar sample where the observational bias does not 
exclude Box 1 of the Hunter-diagram to be populated by rota- 
tionally mixed single stars — which according to Fig.|4]appears 
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indeed possible. While a population of stars which contradict 
the underlying stellar models (presumably binaries) might still 
be present in such a sample, expectedly it would be diluted in a 
dominant population of single stars. 

5. Results: Class 2 stars 

Stars are of Class 2 if the observed surface nitrogen abundance 
is higher than the predicted value of a star rotating with an un- 
protected velocity equal to the v sin(/) of our star at the end of 
the main sequence. This is the case for most of our sample stars 
residing in Box 3 of the Hunter-diagra m (see Fig . [8l>, wh ich ac- 
cording to the statistical analysis of iBrott et al.l (1201 lbl) agree 
reasonably well with the single star evolution models including 
rotational mixing. None of the LMC Class 2 stars is located in 
Box 1 of the Hunter-diagram. While we can not derive any age 
constraints, we can here, relying on single star models with ro- 
tational mixing, constrain their inclination angles. 

The results of the analysis for our sample are summarized 
in Tables IB. II and IB. 21 where we give the upper limits on the 
inclination angle for all stars. Figure ITU1 depicts the histograms 
for the nitrogen inclination angles for the 28 LMC and 24 SMC 
Class 2 stars of our sample. Assuming random orientations of 
stars, the distribution of sin(/) increases for increasing values of 
sin(/). Since we derive upper limits on the inclinations for stars 
with determined surface nitrogen abundance, the distribution of 
those can be expected to be shifted to even higher values. 

The expected behaviour is seen for our LMC stars, where the 
distribution of the sine of the nitrogen inclination angles favours 
high values with a maximum between 0.7 and 0.8. We refrain 
from a further analysis of this histogram since the selection of 
Class 2 stars introduces a bias against the fastest rotators (see 
Fig.|4|. While the number of very fast rotators is small, this in- 
troduces a bias in our distribution which is difficult to account 
for. 

However, it is interesting to note that the distribution of the 
nitrogen inclination angles of the SMC stars shows a different 
behaviour. Here, high values are rare and the average sin(iN) 
value is about 0.5. 

For most SMC stars, only an upper limit to the surface nitro- 
gen abundance is known. The higher his limit, the lower is the 
sin(iN) value required to reach this value. Thus, if the true sur- 
face nitrogen abundance is below the value of the upper limit, the 
sin(iN) value derived from the upper limit will be smaller than 
the one which would be derived from the true nitrogen abun- 
dance. This implies that the sin(/N) value derived from an upper 
limit on nitrogen may be smaller than the true sin /-value of the 
star. Therefore the maximum, or mean, of the sin(iN) distribu- 
tion can be located at lower sin(/)-values than the maximum of 
the true sin(/) distribution. 

Interpreting the distribution of nitrogen inclinations of the 
SMC stars in Fig. [10] in this way, and assuming randomly ori- 
ented rotational axes of these stars, we conclude that the nitrogen 
upper limits derived for many of the SMC stars are significantly 
larger than their true nitrogen abundances. Considering Fig. 7 of 
Hunter et al. (2009), which compares the location of the SMC 
stars with evolutionary tracks in the Hunter-diagram, one can see 
that such a significant shift would indeed be required to bring the 
observations into agreement with the models. 

Similar to the analysis of Class 1 stars, we can use the 
isochrone ages of Class 2 stars with determined nitrogen abun- 
dance to derive their inclination angles. If only an upper limit 
to the surface nitrogen abundance is known, a lower limit to the 
true inclination angle can be found. These inclination angles are 



designated as /isochrone and their sine is given in Tables IB. II and 
IB.2l for all sample stars. 



6. Conclusions 

We present a new method of nitrogen chronology, which can 
constrain the age of a star, and/or its inclination angle, based 
on its observed surface nitrogen abundance, mass and projected 
surface rotational velocity, by comparing the observed nitrogen 
abundance with the one predicted by the theory of rotational 
mixing in single stars. 

This method can be applied to stars when their surface ni- 
trogen abundance increases monotonically with time during the 
main sequence evolution. It has been worked out here for stars 
in the mass range between 5 M and up to 35-50 M Q , and for 
metallicities adequate for the Milky Way the LMC and the S MC, 
based on the stellar evolution models of Br ott et al.l (1201 lab . We 
apply our method to 79 stars from the early B type LMC and 
SMC samples of the VL T-FLAMES Survey of Massive Stars 
(lEvans et al .1120051 l2006h . 

Age constraints from nitrogen chronology could be obtained 
for 17 of the 79 analyzed stars. In Sect. 4.2, we compared those 
to ages of these stars as derived from isochrone fitting in the HR 
diagram (Fig.|7j. We found the isochrone ages of 10 objects to be 
incompatible with their nitrogen age constra ints. Based on their 
rapid rotation and low nitrogen enrichment, iBrott et all (1201 lbl) 
concluded that these 10 stars did not evolve as rotationally mixed 
single stars. Our star-by-star analysis of these objects is reinforc- 
ing their conclusion that the weakly enriched fast rotators in the 
LMC early B stars of the FLAMES Survey — which are 15% of 
the survey stars — can not be explained by rotating single stars. 

For the remaining 7 stars with nitrogen age constraints, six 
LMC and one SMC star, nitrogen and isochrone ages are found 
to be consistent within error estimates. Four of them show good 
agreement, which also allows to derive their inclination angles 
(Sect. 4.2). 

We conclude that our nitrogen chronology results are con- 
sistent with the isochrone method for stars which have no indi- 
cation of an unusual evolution, but incompa tible for stars which 
are suspected binary product according to Brot t et al.l (t201 lbl) . 
Therefore, our new method provides important new results for 
both groups of stars. However, it is clearly desirable to obtain 
larger samples of stars especially of the first group. 

For the 62 stars (28 LMC and 34 SMC stars) for which no 
age constraint could be derived, our method provides limits on 
the inclination for each object (Sect. 5). While the results for 
the LMC stars are roughly consistent with random inclination 
angles, this appears to be different for the SMC stars, for which 
only u pper limits to the n itrogen abundance are available in most 
cases (Hunter et al .120091) . We argue that random inclinations are 
also realized for the SMC stars, with the implication that the true 
nitrogen surface abundances of most SMC stars are significantly 
below the derived upper limits. 

We believe that our new method can be used in the near fu- 
ture to provide further tests of the theory of rotational mixing in 
stars, and — if such tests converge to confirm this theory — to 
act as a new chronometer to constrain the ages of massive main 
sequence stars. 

Based on the method presented here, a web tooQ can be 
found online. For given mass, (projected) surface rotational ve- 
locity and metallicity, the surface nitrogen abundance is calcu- 



1 http://www.astro.uni-bonn.de/stars/resources.html 
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Fig. 10. Histograms of the sin(/Ar)-values for the Class 2 stars of our sample located in the LMC (left panel) and SMC (right panel). 
The number of stars within an area of A sin(iAi) = 0.1 were counted. 



lated and compared to the observed surface nitrogen abundance 
to constrain the age and inclination angle. 
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Appendix A: Equations 

A. 1. Determining parameters a, b, c and d (Eq. @) for the 
example of the LMC metallicity 

In this section, we derive the parameters a, b, c and d of Eq. (f2j). 
This is done for the example of the LMC grid of stellar evolution 
models presented in lBrott et al.l(l201 lal) . 

Our dataset of the stellar evolution models contains informa- 
tion about the time-evolution of the surface nitrogen abundance. 
We can derive the main sequence lifetime from our models and 
therefore obtain the surface nitrogen abundance as a function of 
the fraction of the main sequence lifetime. 

The parameter a corresponds to the initial surface nitrogen 
abundance and is set to be the surface nitrogen abundance at 
t = 0, which for the LMC is given by the constant value 



a = 6.9. 



(A.l) 



By fitting Eq. (0 to the data of the stellar evolution models 
of the LMC grid, we obtain the remaining three parameters for 
different masses and surface rotational velocities at the ZAMS. 

In Fig. IA. H the fitted parameter b is shown as a function of 
vzams for several different initial masses. A polynomial function 
of second order reproduces the data for fixed mass well. The 
parameter b is therefore calculated as 

b = bi ■ vl AM5 + b 2 ■ vzams + b 3 . 



Figure lAT2l shows the coefficients b\, b 2 and bj, as a function 
of initial mass. They can be fitted well by linear approximations 
shown by the solid lines. Using b\(M), b 2 (M) and b 3 (M) the 
parameter b (M [M Q ] , vzams [km/s]) can be calculated using 



M + 4 - 6 f tlji a t\ V ZAMS 

■ +{M + 3.2)- 



2.0 



131 



(M + 0.8).(^) 2 .(A.2) 



Similarly we obtained parameter c as a function of initial 
mass and surface rotational velocity by fitting Eq. (f2} to the data 
of the stellar model grids. The result is shown in Fig. IA.3l on the 
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Fig. A.l. The parameter b is shown as a function of initial sur- 
face rotational velocity for several different masses, obtained by 
fitting Eq. (0 to the stellar evolution grid data (LMC). The solid 
line represents the best fit for the 30 M models. 



left panel. The data is fitted well as a function of Vzams using the 
polynomial function of second order 



C = C\ ■ V 



ZAMS 



+ c 2 ■ Vzams + C3. 



The coefficients c\, c 2 and C3 are obtained for different initial 
masses to consider a mass dependence. 

In the case of LMC models the dataset seen on the left 
panel in Fig. IA.3l can be reproduced well by a linear function. 
Therefore parameter c\(M) is set to zero. This is not the case for 
SMC and MW were the polynomial function need to be used. 

Figure IA.3I shows the mass dependent coefficients c 2 (M) 
and c^(M) which can be fitted well using linear functions 
shown by the solid lines. Using c\{M), c 2 (M) and ci(M), 
c (M [M G ] , vzams [km/s]) can be calculated by 



M-29.5 vzams 
c = — + (M + 7.3) ■ 



54 



10455 



(A.3) 



Parameter d is obtained for different masses and surface ro- 
tational velocities at the ZAMS. The resulting values are shown 
in Fig. IA.4I It can be seen, that only a weak mass dependence 
exists in the case of our data. For simplification, we neglect the 
mass dependence. Therefore a linear approximation is used to 
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Fig.A.2. The three coefficients bi(M), biiM) and b^M) of the polynomial function to describe the parameter b are depicted as 
functions of mass. All parameters can be fitted well using linear functions which are plotted by solid lines. 
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Fig. A.3. On the left panel, the parameter c is shown as a function of initial surface rotational velocity for several different masses. 
The line represents the best fit for the 30 M Q model. On the right panel, the coefficients C2(M) and ct,{M) are plotted as functions of 
the initial mass. They can be reproduced best using linear functions shown by solid lines. 
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Fig. A.4. The parameter d is shown as a function of initial sur- 
face rotational velocity for several different masses, obtained by 
fitting Eq. ([2]) to the stellar evolution grid data (LMC). The line 
represents the best fit using a linear approximation. A small mass 
dependence can be seen, which will be neglected for simplifica- 
tion. 

describe parameter d as depicted in Fig. IA.4I which can be de- 
scribed by 

d = 9l _vzws (A4) 



A.2. The main sequence lifetime and the surface rotational 
velocity 

The fraction of the main sequence lifetime is used in Eq. (|2]i to 
calculate the surface nitrogen abundance. To be able to constrain 
the absolute age of the star it is necessary to know the time a star 
spends on the main sequence for given initial parameters. 

We therefore derive a formula for the main sequence life- 
time for stars of a given initial mass and surface rotational ve- 
locity with MW, LMC and SMC metallicity. Figure IA.5I shows 
the main sequence lifetime as a function of the initial surface ro- 
tational velocity for the example of the LMC data for different 
initial masses. The higher the initial surface rotational velocity 
the higher the main sequence lifetime. Considering Fig. IA.5l it is 
noticeable that the main sequence lifetime is a strong function of 
the initial mass, but shows only a weak dependence on the initial 
surface rotational velocity. The data can be fitted well using 

— =P- (vzams) 9 + r, (A.5) 
yr 

taking mass dependent parameters p, q and r into account. For 
reproducing the data best, q is chosen to be a linear function 
of mass, which in case of LMC and MW can be simplified to a 
constant value, p and r can be calculated depending on the initial 
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Fig. A.5. On the left panel the main sequence lifetime is shown as a function of the initial surface rotational velocity. Stellar evolution 
grid models for LMC metallicity are depicted for several different initial masses from 5 to 45 M Q . The data can be fitted well using 
Eq. ( IA.5b as shown by the additional lines. On the right panel the two coefficients p and r are depicted as a function of mass. The 
dotted lines represent the best fit using Eq. 



mass. We therefore have 

P = Pi- M P2 + p 3 , 
q — qi ■ M + cj3, 
r - n ■ M n + r 3 . 

By fitting the parameters p, q and r, we gain the following 
formula to calculate the main sequence lifetime: 



yr 



'AT 2 +4.9804- 10" 
v 6.4689 • 10 3 
AT 2 + 1.3251 ■ 1(T 3 



' (VZAMS) 



5.1564-10 



-10 



(A.6) 



In our stellar evolution m odels, mass loss is included as de- 
scribed in IVinket al.l(l2010h . During the main sequence evolu- 
tion the stars loses angular momentum due to mass loss. The 
value of the surface rotational velocity therefore changes as a 
function of time. Figure lA~6l depicts the surface rotational veloc- 
ity of a few exemplary models as functions of time (left panel) 
and effective temperature T e s (right panel). An initial surface ro- 
tational velocity of about 170 km/s is chosen. 

During the evolution on the main sequence the surface ro- 
tational velocity of a star remains almost constant. At the end a 
sudden drop is visible for the highest considered masses, caused 
by an increase in mass loss and therefore a decrease of angular 
momentum when the effective tem perature falls below 25 000 K, 
caused by the bi-stability braking (IVink et al.ll2000ll200lh . Stars 
with masses up to 30 M Q reach the end of the main sequence at 
temperatures higher than 20 000 K (see Fig. IA.6b . They do not 
stay long enough on the main sequence at T e s < 25000 K to 
be significantly slowed down. In the case of stars more massive 
than 30 M the assumption is valid for T e ff > 25000 K. For stars 
which have not been influenced noticeably by this effect it is rea- 
sonable to estimate the surface rotational velocity by the value 
at the ZAMS. 



A.3. Small Magellanic Cloud 

For the SMC, the parameters specifying Eq. (|2) are: 

a = 6.5, 



M + °- 9 + (M + 3.4), VzAMS 



1.4 

M-71 



122 



<« + 2.0).(l^) 2 



c = 

d = 9.0 



yr 



67 

_ I'ZAMS 

~ 58 



+ 40»-^ + (M + 213)-f^) , 



18457 



V 6294 



M 



-5.85 



(V 



ZAMS 



0.1121-M+2.4187 



+ - 



8.5881 ■ 10- 4 

M -2.065 + L6507 . 10 -3 



5.0320 • 10 



■10 



The formula to calculate the main sequence lifetime for SMC 
stars is valid for all stars analyzed, but not in the complete range 
given in Table [2] Table IA.1I summarizes the mass and velocity 
ranges valid for this specific formula. 



Table A.l. Range of initial stellar parameters for which the main 
sequence lifetime calculation for SMC stars is validated. 
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A.4. Milky Way 

For the MW, the parameters specifying Eq. (f2) are: 
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Fig. A.6. The evolution of the surface rotational velocity over time is shown on the left panel, while the surface rotational velocity as 
a function of the effective temperature is given on the right panel. Models of the stellar evolution grid in the case of LMC metallicity 
are used for different initial masses and surface rotational velocities of about 170 km/s. 
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Table B.l. Properties of the analyzed LMC stars. 



Star 


Number 


v sin (0 




M cvo l 


N 




sin (;'n) 


sin (iiso) 


f]M 


'isochrone 


binary 


class 






km/s 


K 




dex 


Myr 






Myr 


Myr 






N 1 1-34 


LI 


203 


25500 


20 


7.12 


7.5 


1.00 


- 


^•'-1.5 


7.8 


yes 


1 


N 1 1-37 


L2 


100 


28100 


23 


<7.17 


6.2 


0.76 


>0.76 


6.2 


6.8 


yes 


2 


N 11^6 


L3 


205 


33500 


28 


<7.71 


5.1 


0.88 


>0.88 


5.1-2.5 


4.8 


yes 


2 


N 1 1-77 


L4 


117 


21500 


12 


6.96 


16.0 


0.90 


0.89 


16 ifin 


15.4 


yes 


2 


N 11-88 


L5 


240 


24150 


14 


<6.96 


12.6 


1.00 




2 6 +1 - 9 


12.6 




1 


N 1 1-89 


L6 


117 


21700 


12 


6.91 


16.0 


1.00 


- 


14 6 +1 4 


16.4 


yes 


1 


N 1 1-95 


L7 


267 


26800 


15 


7.46 


11.4 


1.00 


- 


4 0+2-* 


10.4 


- 


1 


N 11-102 


L8 


218 


31000 


18 


<7.58 


8.6 


1.00 




6.6+|° 


7.0 


— 


1 


N 11-104 


L9 


153 


25700 


14 


<7.20 


12.5 


0.99 


>0.97 


12.5-43 


11.8 




2 


N 11-111 


L10 


101 


22950 


12 


6.93 


16.0 


0.83 


0.81 


16.0_ 1M 


15.4 


- 


2 


N 11-114 


Lll 


299 


32500 


19 


<7.92 


8.2 


0.91 


>0.90 


8-2-4.7 


5.8 


- 


2 


N 11-116 


L12 


160 


21700 


11 


6.98 


18.6 


1.00 


- 


12 4 +62 


18.4 


- 


1 


N 11-117 


L13 


247 


26800 


13 


<7.06 


14.2 


1.00 




4 r 1 - 9 

^• x -4.1 


11.2 




1 


N 11-118 


L14 


150 


25700 


13 


7.39 


14.1 


0.83 


0.79 


14.1 


12.6 


yes 


2 


N 11-120 


L15 


207 


31500 


17 


<7.56 


9.3 


1.00 


0.91 




6.4 




1 


N 11-121 


L16 


265 


27000 


14 


<7.48 


12.6 


1.00 


- 


6-1^ 


11.2 


- 


1 


N 11-122 


L17 


173 


33000 


18 


<7.72 


8.6 


0.68 


>0.61 


8.6 


4.6 


- 


2 


NGC 2004-20 


L18 


145 


22950 


18 


7.46 


8.6 


0.83 


0.83 


8.6 


9.6 


yes 


2 


NGC ?f)04-30 


L19 


123 


29000 


19 


7.43 


7.9 


0.71 


0.71 


7.9 


7.6 




2 


NGC 2004-47 


L20 


133 


21700 


12 


7.03 


16.0 


0.94 


0.94 


16.0-i 6 .o 


16.4 


yes 


2 


NGC 2004-50 


L21 


109 


20450 


11 


7.10 


18.6 


0.71 


0.68 


18.6-18.6 


17.0 


yes 


2 


NGC 2004-54 


L22 


114 


21700 


12 


6.97 


16.0 


0.86 


0.86 


16.0_ 160 


16.4 


yes 


2 


NGC 2004-62 


L23 


106 


31750 


17 


<7.33 


9.3 


0.67 


>0.61 


9.3 


7.6 


- 


2 


NGC 2004-63 


L24 


107 


21700 


11 


6.97 


18.6 


0.79 


0.79 


18.6-18.6 


18.0 




2 


NGC 2004-65 


L25 


223 


20450 


11 


<7.19 


18.7 


1.00 




9 7 +3 - 7 

y - '-3.6 


19.8 




1 


NGC 2004-66 


L26 


238 


25700 


13 


7.98 


14.2 


0.60 


0.59 


14.2 


12.0 




2 


NGC 2004-69 


L27 


178 


28000 


15 


7.84 


11.2 


0.57 


0.56 


11.2 


9.8 




2 


NGC 2004-74 


L28 


130 


27375 


14 


<7.42 


12.5 


0.72 


>0.67 


12.5 


10.6 


yes 


2 


NGC 2004-75 


L29 


116 


21700 


11 


7.06 


18.6 


0.78 


0.77 


18.4_,8.6 


18.0 




2 


NGC 2004-77 


L30 


215 


29500 


16 


7.65 


10.2 


0.93 


0.90 


10.2-3.6 


8.4 




2 



NGC 2004-79 
NGC 2004-81 
NGC 2004-82 
NGC 2004-85 
NGC 2004-88 


L31 
L32 
L33 
L34 
L35 


165 
105 
161 
150 
200 


21700 
26800 
25700 
20450 
20450 


11 

13 
13 
10 
10 


7.60 
7.14 
7.47 
7.25 
<7.18 


18.6 
14.0 
14.1 
22.0 
22.0 


0.70 
0.70 
0.83 
0.85 
1.00 


0.70 
0.63 
0.79 
0.84 


18.6 
14.0 
14.1 
22.0 

15 r" 


18.4 
11.2 
12.6 
21.4 
21.6 


yes 
yes 


2 
2 
2 
2 
1 


NGC 2004-95 
NGC 2004-99 
NGC 2004-100 
NGC 2004-101 
NGC 2004-104 


L36 
L37 
L38 
L39 
L40 


138 
119 
323 
131 
274 


25700 
21700 
26800 
21700 
25700 


13 
10 
13 
10 
12 


7.13 
7.06 

<7.29 
7.59 

<7.28 


14.1 

22.0 
14.5 
22.0 
16.3 


0.92 
0.78 
1.00 
0.55 
1.00 


0.89 
0.74 

0.53 


14.1-4.8 

22.0_ 22 .o 
3 6 +L3 
22.0 


12.8 
19.6 
11.6 
19.6 
13.2 




2 
2 
1 

2 
1 


NGC 2004-105 
NGC 2004-107 
NGC 2004-110 
NGC 2004-113 


L41 
L42 
L43 
L44 


235 
146 
121 
307 


25700 
28500 
21700 
20450 


12 
14 
10 

9 


<7.11 
<7.12 
6.95 
<7.29 


16.2 
12.5 
22.0 
27.1 


1.00 
1.00 
0.91 
1.00 


0.88 
0.87 


6.1!" 
12 3 +0 - 2 

22.0-22.0 

9 5 +36 


13.2 
9.4 
19.4 
22.2 


yes 


1 
1 

2 
1 



Notes. From left to right: catalogue number (Evans et al. 2005), reference number used in this paper, projected surface rotational velocity, effective 
temperature, evolutionary mass and surface nitrogen abundance (all from Hunter et al. 2008b, 2009), main sequence lifetime, nitrogen inclination 
angle, inclination angle derived using the nitrogen abundance and the isochrone age, maximum age based on the surface nitrogen abundance, 
isochrone age, identification of binarity and class. For Class 2 stars, adopting an error in the nitrogen abundance does not lead to an increased 
upper age limit, and in many cases also not to an decreased lower age limit. 
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Table B.2. Properties of the analyzed SMC stars. 



Star 


Number 


v sin (0 


Teff 


M cvol 


N 




sin 0' N ) 


sin 0' iso ) 


tN 


^isochrone 


binary 


class 






km/s 


K 


M e 


dex 


Myr 






Myr 


Myr 






NGC 346-27 


SI 


220 


31000 


18 


<7.71 


9.7 


0.57 


>0.57 


9.7 


7.6 


- 


2 


NGC 346-32 


S2 


125 


29000 


17 


<6.88 


9.1 


0.85 


>0.85 


9.1 


9.2 


yes 


2 


NGC 346-53 


S3 


170 


29500 


15 


<7.44 


10.9 


0.57 


>0.57 


10.9 


9.4 


yes 


2 










i j 


<7 9S 


1U.7 




^_U. JO 


1 n o 

1U.7 


Q 4 




9 


NGC 346-58 


S5 


180 


29500 


14 


<7.56 


12.0 


0.51 


>0.51 


12.0 


9.6 


yes 


2 


NGC 346-70 


S6 


109 


30500 


15 


<7.50 


11.0 


0.34 


>0.33 


11.0 


8.2 


- 


2 


NGC 346-79 


S7 


293 


29500 


14 


<7.88 


13.4 


0.65 


>0.64 


13.4 


9.4 


- 


2 


NGC 346-80 


S8 


216 


27300 


12 


<7.85 


15.4 


0.48 


>0.48 


15.4 


12.8 


- 


2 






Z.J J 


Z. 1 zuu 


q 


<7 ^9 


9S 1 


n 73 


>n 73 


9^ 1 

Z J . 1 


94 4 






NGC 346-82 


S10 


168 


21200 


9 


<7.55 


24.7 


0.46 


>0.46 
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Notes. From left to right: catalogue number (Evans et al. 2005), reference number used in this paper, projected surface rotational velocity, effective 
temperature, evolutionary mass and surface nitrogen abundance (all from Hunter et al. 2008b, 2009), main sequence lifetime, nitrogen inclination 
angle, inclination angle derived using the nitrogen abundance and the isochrone age, maximum age based on the surface nitrogen abundance, 
isochrone age, identification of binarity and class. For Class 2 stars, adopting an error in the nitrogen abundance does not lead to an increased 
upper age limit, and in many cases also not to an decreased lower age limit. 
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